Splitting of Landau levels of a 2D electron due to electron-phonon interactions 
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We show that in a very strong magnetic field B electron-phonon interaction gives rise to a splitting 
of Landau levels of a 2D electron into a series of infinitely degenerate sublevels. We provide both 
qualitative and quantitative description of this phenomenon. The cases of interaction with acoustic 
and polar optical phonons are considered. The energy distance between nearest sublevels in both 
cases tends to zero as B~ 1/2 at large B. 



The dynamics of a 2D lowest Landau level (LL) elec- 
tron in high mobility samples where only the electron- 
phonon interaction is relevant was disscussed in Refs. 
jl],^| . In the limit of weak electron-phonon coupling the 
temperature dependence of the longitudinal electron mo- 
bility was calculated in these works both at low |jj and 
relatively high 0] temperatures. 

The limiting case of weak electron-phonon interaction 
is not, however, realized at very high magnetic fields. The 
effective coupling constant a characterizing the electron- 
phonon interaction in a strong magnetic field increases 
with B (see below) and may reach large values at which 
a new physics comes into play due to the essential role 
of polaronic effects. In this limit the conclusions of Refs. 

are no longer valid. Some aspects of the problem 
of calculating the electron mobility in this case were dis- 
cussed in Ref . || . 

In this letter we would like to concentrate on a simpler 
problem of finding the energy spectrum of the electron- 
phonon system when B — > oo. Knowledge of this spec- 
trum is nesesary for the correct calculation of the electron 
mobility at very low temperatures. 

For ci» 1 the lattice in the vicinity of an electron is 
strongly deformed and this deformation creates a pola- 
ronic well (PW) around the electron. The characteristic 
frequency of the electron motion in the PW is very high 
compared to the phonon frequencies and the lattice de- 
formation forming the PW can be considered quasistatic. 
The energy and the wave function of the ground state 
(GS) of the electron in the PW in the limit B — ► oo were 
calculated in Ref. jlj in the framework of Pekar adiabatic 
method JH. At first sight it may seem that the results 
of Ref. H| give the complete description of the physics 
of the problem and the polaronic effects simply reduce 
to an energy shift of the LL. This, in fact, is not true. 
As we will show, in strong coupling limit the LL splits 
into a series of sublevels, these sublevels being infinitely 
degenerate like the original LL. 

This splitting can be understood qualitatively on an 
example of a simple model, which can be solved exactly. 
Let us consider a 2D electron with the effective mass 
to* placed in a strong magnetic field and coupled by a 
harmonic potential to a heavy neutral particle with the 



mass M 3> to*. This particle models the phonon cloud 
surrounding an electron in the crystall. There are three 
different types of motion in the system. The fastest one is 
the motion of the electron with the cyclotronic frequency 
lu c = eB/m*c. The second one is a slower drift of the 
magnetic circle center of the electron around the heavy 
particle with the frequency fl = u>q/u) c , where loq <C to c is 
the frequency of the electron oscillations in the harmonic 
potential in the absence of the magnetic field. The third, 
slowest type of motion is the cyclotronic rotation of the 
system as a whole with the frequency fi c = eB/Mc. The 
energy spectrum of the model can be classified by three 
quantum numbers n, to, N corresponding to these types 
of motion E n ^ m ^ = fko c n + Mlm + hQ c N . 

In case of an electron interacting with phonon field 
there also exist three types of motion. In analogy with 
the model, the Landau levels corresponding to the fast 
cyclotronic rotation of the electron are split due to the 
slower drift of the magnetic circle center in the PW and 
to yet slower rotation of the polaron as a whole. In this 
letter we restrict ourselves to the calculation of the split- 
ting of the lowest LL. As we will show the drift frequency 

depends on the radius of the drift orbit and the series 
associated with the quantum number to is not equidis- 
tant. The frequency of the slowest cyclotronic rotation 
of the polaron as a whole Q c is inversely proportional to 
the polaronic mass M, which depends on the magnetic 
field and the quantum number to. The appearance of 
a finite mass of the polaron (an electron at the LL has 
effectively no mass) , is a result of dressing of the electron 
by phonons. We want to emphasize that in our consider- 
ations the transitions between different LLs of the elec- 
tron are neglected and the finite polaronic mass arises 
when the states of one LL only are taken into account. 
In strong coupling limit the mass M is proportional to a 
and grows with magnetic field as B 3 / 2 , whereas the dis- 
tance between nearast sublevels Ml c decreases as Z? -1 ' 2 . 
Note also, that all the sublevels of the spectrum obtained 
in our letter are infinitely degenerate over the position of 
the magnetic circle center of the polaron as a whole. 

We consider an electron at zero temperature in a quan- 
tum well (QW) of the width a located in the plane (x, y) 
to which a strong magnetic field B parallel to the axis 



1 



z is applied. We assume that the electron is at the first 
level in QW and at the lowest LL. The sample is sup- 
posed to be fairly pure and the filling factor of the LL is 
supposed to be small so that the impurity scattering as 
well as the electron-electron scattering can be neglected. 
The Hamiltonian of the system projected on the low- 
est LL can be written |^|| in terms of the operator of 
magnetic circle center R = (X,Y) with non-commuting 
components [X, Y) = il 2 : 



H 
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/ q Q q sin(q||R-7r/4), (1) 



where I is the magnetic length, Q q are the normal coordi- 
nates of the oscillations of the lattice and P q are the mo- 
menta conjugated to them, w q is the phonon frequency, 
q is the phonon wave vector with the components q| | in 
the QW's plane and q z in z direction. The function / q 
can be written as 



V (q z ) = / dze l ** z \Mz)\ 2 - (2) 



/ q = FqTj(q z )e~ 



The factors exp(— g 2 Z 2 /4) and f](q z ) appear as a result 
of the projecting of the original Hamiltonian of the sys- 
tem on the zero LL and on the lowest level in the QW 
with the wave function of transversal quantization ^po(z). 
Due to these factors the electron effectively interacts 
with the phonon modes with momenta lying in the shell 
<7|| < l/l, q z < 1/a. The function F q is proportional 
to the square root of coupling constant a. Its concrete 
form depends on the type of phonons. We will do all the 
calculations for arbitrary and afterwards consider two 
examples of polar optical phonons and acoustic phonons 
interacting with the electron via deformation potential. 

We assume that at a > 1 characterisic frequency of 
electron motion in the PW is much larger than phonon 
frequencies. Then the PW is quasistatic and one can in- 
troduce the wave function ^ of the GS of the electron 
in this well. Averaging the Hamiltonian ([!]) with W we 
get the Hamiltonian which only depends on the phonon 
coordinates and momenta. The ground state energy of 
this Hamiltonian is given by 



f 2 

q 1 



= - ^ ^<*| sin(q N R - n/4)\V) 2 . 



(3) 



In Pekar method |!| \& should be obtained by minimizing 
the energy functional (||) under the constraint = 1. 

One can see that the functional (||) is stationary on the 
eigenf unctions of the angular momentum operator of the 
electron at the lowest LL M = h[l 2 - X 2 - Y 2 ]/2l 2 . 
This implies that the self consistent potential of the 
PW is axially symmetric. The absolute minimum of 
the energy functional (|3|) is reached on the zero angu- 
lar momentum wave function [Q which in X represen- 
tation ( X = X and Y = —U 2 d/dX) has the form 



V (X) = (7rP)- 1 /4 exp (_x 2 /2/ 2 ). The GS energy is 
given by 



(4) 



Now let us consider the excited polaronic states in 
the quasiclassical approximation. The Lagrange func- 
tion corresponding to the Hamiltonian ([!]) is written as 
follows 



E^- 



q 



/q<2 q sin(q||R- tt/4). 



(5) 



The equations of motion given by this Lagrange function 
read 



R 



E/q^qt 11 x cos(q||R- tt/4), 



Qq+^q<3q= ~/q Sill(q| |R - 7r/4), 



(0) 
(7) 



where n is the unit vector in z direction. Classical mo- 
tion of the electron in the axially symmetric PW is simply 
rotation 



R — Rn 



(.Rocos(fit) +x, i?osm(Oi) +y), (8) 



where r is the position of the PW. After substituting (@) 
into (0) we see that the phonon coordinates Qq can be 
represented as a sum of harmonics with the frequencies 
which are integer multiples of Q. Assuming that fi ^> cj q 
we can keep in the sum only the harmonic with zero fre- 
quency and obtain the following expression for Q q 



Qq * Q>) 



Jo(q\\Ro)fq . , , 
H sin(qi | r - tt/4 . (9) 



Substituting Eqs. @ and (||) into @ we get the fre- 
quency f2 of the classical drift of the electron in the ef- 
fective phonon potential U(Rq) 



n 



l 2 dU{R Q ) 

HRq C?i?o 



U(Ro) 



q 



(10) 



Quasiclassically the values of the rotation radius Rq are 
quantized as 



Ran 



l\/2m, m ^> 1, 



(11) 



where m is an absolute value of the quantum of the angu- 
lar momentum. In our case m coinsides with the number 
of the energy level in the PW. The energy of the m-th 
high excited state is given by 



U(R t 



o in } 



(121 



These levels are separated by large energy ~ h£l propor- 
tional to the coupling constant. 



2 



The high frequency motion described above is con- 
nected with the internal degrees of freedom of the po- 
laron. Besides these internal degrees of freedom there 
exist two more degrees of freedom x, y (see Eq.(ph) corre- 
sponding to the motion of the polaron as a whole. Now 
we will show that polaron aquires a finite mass M ~ a 
and moves as a heavy particle with the electric charge e 
in magnetic field B. Quantization of this motion leads 
to additional equidistant splitting of energy levels. The 
energy of splitting is equal to the cyclotronic energy 
hfl c = heB/Mc. 

First we consider the problem quasiclassically, assum- 
ing that the electron in the PW is in one of the highly 
exited states ( m ^> 1) . Let us treat vector r in (g) as 
a slowly changing function of time. Then the deforma- 
tion field Q q will adiabatically follow the motion of the 



polaron according to Eq. 
change of variables Q q — 



Q° q (r) 



Let us do the following 
9 q , where 9 q stands 



for the oscillations of the phonon coordinates around Q° . 
Substituting Eqs.(||) into (||) and averaging over the fast 
rotation with the ferquency f2 we get the effective La- 
grange function corresponding to the slow motion of the 
polaron 



£ e £ = -p yx 



M m r 2 



E 



9 q - W q#q 



(13) 



The first two terms in Lagrange function (|13|) describe 
the motion of a massive particle with the electric charge 
e in the magnetic field B. Classically this motion is cy- 
clotronic rotation with the frequency Q c . The mass M m 
of the polaron in m-th exited state is given by 



M„ 



/q^o(<?||-Rom) 



q 



(14) 



The term £; nt in the Lagrange function describes the 
interaction of polaron with the phonons # q . This inter- 
action leads to decay of the cyclotronic rotation of the 
polaron and Q c acquires imaginary part T. For acous- 
tic phonons T is small as compared to O c if the velocity 
of the cyclotronic motion is less than the sound veloc- 
ity. For optical phonons the decay appears as long as 
we take into account the dispersion of phonon frequency. 
For weak dispersion T is also small. In both cases the 
quantization of cyclotronic rotation of the polaron leads 
to well defined Landau levels of the polaron as a whole. 

Quasiclassical approximation is not valid if we con- 
sider the cyclotronic splitting of the GS. Still in this 
case it is also possible to derive the effective Lagrange 
function for the slow motion of the polaron. This can 
be done using the standard procedure of integrating out 
fast components Ro(i) of the vector R(t) in functional 
integral describing transition amplitudes of the electron 
interacting with the phonon field. The resulting effective 
Lagrange function will only differ from ([l3]) by concrete 
form of the interaction term and by the expression for 



the mass M. Namely, in Eq. (|14|) the Bessel function 
Jo(<Z||i?om), which had appeared as a result of classical 
averaging over fast rotation should be replaced by the 
quantum average of exp(«q||R) over the GS of the elec- 
tron exp(iyX/l 2 )~$>o(X — x) in the PW placed at the 
point (x,y). For the mass of the polaron in the GS we 
have 



Mn = 



q 



exp 




(15) 



The highly exited states of an electron in the PW can 
decay into lower exited states, the energy difference being 
transfered to the polaron as a whole. The rigorous cal- 
culation of corresponding transition probability is cum- 
bersome. One can argue, however, that this probability 
is negligibly small. Indeed, the momentum P transfered 
to polaron in this process is P = y/lM(E m — E m i) ~ 
V MhQ, ~ a. In strong coupling limit P is large compared 
to h/l and cannot be compensated byphonon emission 
due to the factor exp(— g^Z 2 /4) in Eq.(g). 

The equations obtained above are general and can be 
applied to any type of electron-phonon interaction. Next 
we consider two concrete examples of an electron inter- 
acting with acoustic and polar optical phonons. For cer- 
tainity we will assume that a <C I. 

1. Acoustic phonons. For acoustic phonons 



Fa 



C' 2 



V 



npM 2 s 3 



(16) 



Here a ~ B is the dimcnsionless coupling constant char- 
acterizing the electron-phonon interaction in a strong 
magnetic field , C is the deformation potential constant, 
p is the density of the crystal, s is the sound velosity, and 
V is the volume of the crystall. 

Performing the summation in Eqs.(Q) we obtain energy 
of the GS of the polaron 



En = — 



afioj a 



s 



dq zV 2 (q z ) ~ - (17) 
a 



Classical frequency of rotation of the electron in the 
PW is determined by Eq.([Io| ). Taking into account that 
m Eq.flfiJ) Ro > I we obtain the following expression for 
the drift frequency 



n 



(18) 



Eq. (|T^) was derived under the assumption that f2 3> u) q . 
The phonon frequencies are restricted by u> a . Taking into 
account Eq. (|lTJ) we have the following condition of the 
validity of our calculations 



a>l, m <C a 



2/3 



(19) 



By the virtue of @ and @ for the energy levels in 
quasiclassical approximation we have 
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Brn, 



ahuj a 
8tt 3 / 2 



(20) 



The mass of the polaron in ground and highly excited 
states are obtained from Eqs. (|14] . [l5| ) 



M = a 



3277' 



M m = a- 



(21) 



The spacing between Landau levels of the polaron de- 
pends on m and is equal to Ml a c — KeB/(Moc) and 
ftfi™ = %eB I (M m c) for ground and highly exited states 
correspondingly. Viewing the term C m t in Eq.(p^[) as 
a perturbation one can get the following result for the 
imaginary parts r and r m of the cyclotronic frequen- 
cies Q° and fit? 



The above results are in contrast with the usual point 
of view that the motion of the electron strongly inter- 
acting with phonons in a magnetic held is characterized 
by two frequencies only: the frequency of the fast motion 
of the electron in the PW and the frequency of the cy- 
clotronic rotation of the polaron as a whole with the mass 
not depending (or slightly depending) on B. In fact this 
is true only when oj c is small compared with the charac- 
teristic frequency of the electron motion in the PW. Our 
results are valid in the opposite case. 

Note finally, that the splitting of the LL can be ap- 
parently observed in Ge samples, where the deformation 
potential is considerably large. Simple estimation shows 
that in this material effective coupling constant a be- 
comes sufficiently large at B > 50 T. 
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(22) 



to the first order in small parameter 1/a 3 . We see that 
in a strong magnetic field when a is large T is small com- 
pared to the spacing f2 c and Landau levels of the polaron 
are well defined. Eqs.(p2|) are valid if the polaron velos- 
ity is smaller than the sound velosity . In the opposite 
case r becomes very large and cyclotronic quantization 
of polaron motion is absent. 

2. Optical phonons. For optical phonons 



F q = 



1 8alhuiQ 



flLOnl 



1 



(23) 
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Here a ~ \J B is the dimensionless coupling constant for 
optical phonons in a strong magnetic field and ujq is the 
optical phonon frequency . The expressions for the drift 
frequency and the energy of an electron and the mass of 
the polaron in ground and highly exited states read 



Q = 



au)o Z 3 ln(i? /0 



2vr 2 



R 3 



M n = 



ahujQ 



ah 



ahajQ him 



r, M m = 



It follows from Eqs.( pT| , p6[ ) that for given m 
M~B 3/2 , tt r ~B- 1/2 . 
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